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Abstract. In the presence of viscosity the hydraulic jump in one dimension is seen to be a first-order 
transition. A scaling relation for the position of the jump has been determined by applying an averaging 
technique on the stationary hydrodynamic equations. This gives a linear height profile before the jump, as 
well as a clear dependence of the magnitude of the jump on the outer boundary condition. The importance 
of viscosity in the jump formation has been convincingly established, and its physical basis has been 
understood by a time-dependent analysis of the fiow equations. In doing so, a very close correspondence 
has been revealed between a perturbation equation for the fiow rate and the metric of an acoustic white 
hole. We finally provide experimental support for our heuristically developed theory. 

PACS. 47.15. Cb Laminar boundary layers - 47.60.+i Flows in channels - 47.32.Ff Separated flows 



J> 1 Introduction hydraulic jump PI^IE]. It is a familiar observation, seen 

everyday in the kitchen sink. A similarly abrupt increase 



■ A stream of water impinging vertically on to a horizon- 



in the height — a jump — occurs in the one-dimensional 
flow as well, and this phenomenon finds mention in many 
introductory text books on hydrodynamics QlEliniEllHl ■ A 

very regularly cited practical example of the jump in the 

one-dimensional flow is the passage of a tidal bore up a 

river |3. However, the texts include viscosity — arguably 
sbsinghaOtheory . tif r . res . in 

tpjkbSmahendra. iacs .res . in the primary physical cause of the jump — only through a 

arnabOmri . ernet . in 



tal surface spreads out radially in a thin sheet along the 
plane from the point of impingement, and at a certain 
radius the height of the flowing layer of water suddenly 
increases. In this two-dimensional flow, such an abrupt 
increase in the level of the liquid is known as the circular 
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phenomenologically added friction term 6 . Consequently, put forward a heuristic picture of the role of viscosity and 

starting from the Navier-Stokes equation, it has not been find that in its presence the jump attains a finite value at 

possible to predict the position of the jump in terms of the critical Froude number. This is exactly what happens 

the volumetric flow rate. For the two-dimensional flow, on in a first-order transition. In our calculations we also es- 

thc other hand, the role of viscosity has been very clearly tablish a connection between the magnitude of the jump 

taken into account in the works of Bohr et al. |5lllUllllj . and the outer boundary condition. We use steady hydro- 

This has led to a scaling dependence for the position of dynamic equations to obtain the position and the height 

the jump on the volumetric flow rate. The two-dimensional of the jump. In the process we reinforce our conclusion 

problem, however, is sufficiently complicated, and in pre- that the jump is of the nature of a first-order transition, 
dieting the position of the jump in this case, knowledge of 
asymptotic solutions has been necessary to a fair extent. 



Furthermore, we carry out a time-dependent linearised 
perturbative analysis about the steady flow solution. We 
perturb the steady and constant flow rate, and see that 
viscous dissipation, of course, drives the system back to- 
wards stability, but what we have made a note of with 

work. We derive a profile of the height of the liquid layer in , , • , , , . r . ^ 

much greater mterest is that the equation tor perturba- 



Motivated by the methods applied to study the two- 
dimensional flow, and by the results obtained thereof, we 
make a similar analysis of the one-dimensional flow in this 



tion in the flow has a remarkable degree of closeness with 
a metric that implies an acoustic white hole. This anal- 
ogy, coupled with some characteristic time scales obtained 
through the perturbative analysis, enables us to argue for 
a physical basis behind the formation of the jump. 



a one-dimensional open-channel flow by making transpar- 
ent approximations about the nature of the flow. We flnd 
that the position of the jump is sensitive to the vertical 
profile of the velocity field, and we make use of this de- 
pendence to conclude that the profile is far from parabolic 
and resembles more closely a turbulent profile. There ex- 
ists no definite result for the magnitude of the jump. Text It is in a largely heuristic spirit that we have made 
books [oIIdIItIIo] analyse the problem in one dimension, our theoretical foray into the channel fiow problem. To 
from which it can be shown that for hi and h2 being the bolster our arguments, we have therefore brought forth 
heights of the liquid layer before and after the jump re- some experimental evidence in support of our theory. We 
spectively, the ratio h2/hi is unity for the critical value have measured the magnitude of the jump and the height 
of the Froude number J- (i.e for J- — 1), and increases profile of the fiow before the jump. On both counts we find 
smoothly as the Froude number is increased from unity, good agreement with our theory. From our experimental 
This indicates that the jump is a second-order transition, data we can also easily infer that the parabolic profile in 
Contradicting this viewpoint, in our present analysis we the vertical direction has no validity. 
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Fig. 1. The geometry of the experimental set up. The hori- 
zontal flow starts from the edge AB. 

2 Role of Viscosity : A Heuristic Study 

The flow occurs in a channel of width L, with L being 
very much greater than the depth of the liquid layer, as 
has been schematically shown in Fig^ At some point the 
depth changes from hi to h2. We work with a control 
volume which extends from a point before the jump to a 
point after the occurrence of the jump. Ignoring the effect 
of viscosity, the continuity equation can be written down 
as 

Lhiui = Lh2U2 = Q (1) 

in which ui and U2 are the flow velocities before and af- 
ter the jump respectively, while Q is the volumetric flow 
rate. The momentum change pQ {u2 — ui) per unit time is 
brought about by the force due to the pressure difference. 
This gives the balancing condition 

^pgL{hl~hl) ^ pQiu2~ui) (2) 

First writing Ti = h2/hi and the Froude number J- as 
!F — ul/ghi, we combine Eqs.lQJ and ||2J), to get 

n{l+n)-2T=0 (3) 

which can be solved to get 27i = — 1 + Vl + 8J-. If we now 
write T = 1 + with the condition < i9 ^ 1, then to 
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Fig. 2. Control volume of the flow, demarcated by the dotted 
lines. The boundary layer has an average thickness of S. 

first order in i) we will get H = l + 2i3/3, which establishes 
the standard text book interpretation [SUnitZI of the jump 
being a continuous transition as a function of J^. 

We need to modify the above picture in the presence 
of viscosity. The most important contribution of viscos- 
ity will be the formation of a boundary layer. Practically 
speaking, the thickness of the boundary layer increases 
as the fiow progresses along the plane, but in our control 
volume we have constrained the flow to have an average 
thickness of 6, within which the velocity increases from 
zero to ui, while over the depth of hi — S, the flow ve- 
locity remains at a constant value of ui. This state of 
affairs has been schematically represented in Fig|21 After 
the jump, beyond a mixing zone characterized by vortices, 
the flow has a mean speed U2 and an increased depth h2. 
However, it must be noted here that following the jump, 
the flow has been known to be turbulent — something 
that may very readily be appreciated from the analogous 
case of the circular hydraulic jump jl2| — and it becomes 
too much complicated to be thought of simply in terms 
of a boundary layer. Therefore we keep our analysis of 
the control volume confined to the flow region before the 
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jump. Within the boundary layer we make use of an arbi- that reads as 

trary profile i.e. u{z) = uiip{zl6) with (p(^) < 1 for all f , / 1 \ 8 

2t{ — -1\-2T- 

and with ^ itself constrained by the range < ^ < 1. The 



n hi 



2(l-/i)^+/2-l 



continuity equation now reads hi rl 

Lui {hi-6) + L "^'^(i) ^ Lu2h2 = Q (4) from which we can easily see that if <5 = 0, i.e. if the 

effect of viscosity is neglected, then the result given by 

which, following some manipulations, can be rendered as 
Lui [hi -S{1- h)] = Lu2h2 = Q (5) 



Eq.lPJ will be recovered. In this inviscid limit, prescribing 
T = I + leads to H = 1 + e with e = 2i9/3. In the 
where h = /g^ ip{^) d^. In the presence of viscosity, this is presence of viscosity, we once again seek a jump solution 



the modified form of Eq.ljTJ. by writing H = 1 + e with e > for = 1 + t?. In the 

The force balance equation requires a similar modifi- limit ■& — > 0, we will then have the cubic equation 



cation in the momentum flow rate. Before the jump the x 

£3 + 3^2 + 2— (1 - I2) 

amount of momentum entering the control volume per 



^hi 



(I + /2-2/1) 



2) 



5^ 



• • • ' 2f^(l-/i)2 = (10) 

unit time is hf 

Lp (hi — 6) + Lp [ u^tf^ ( — \ dz clear that there can be no positive root of e which is 

^ Jo ^ 

greater than {1 + 12 — 2/i)/(l — I2) since in that case all 



= pQui — LpuiS (Ii — I2) (6) 

where I2 — Jq (p^{S,)d^. The amount of momentum leav- 
ing the control volume per unit time is pQu2 and the force 
due to pressure on the control volume is pLg {h\ — /2 
acting to the right. The force balance equation now be- 
comes. 



terms on the left hand side of Ea. (|10|l will be positive. For 

a Couette profile Ii = 1/2 and I2 = 1/3. This gives e — 

0.5 as an upper limit. For the parabolic profile Ii =2/3 

and I2 = 8/15, giving e ~ 0.43 as a maximum value. On 

the other hand, a continuous transition would imply that 

e = 0. If we use this value of e in Eg. 1101) . we will get the 

pQ {u2 - ui) + LpSul {Ii - I2) = ^pgL {hj - hi) (7) condition 6/hi = (1 + /2 - 2/i)/(l - hf, which, for both 

TT • ii j-i- J ■ irr the Couette and the parabolic profiles should give a value 

Using the continuity condition as expressed m Eq.ljSj, we t- t- b 

of S/hi to be greater than unity. This is physically an 
untenable result, because 6 is the average thickness of the 



il[hi - <5(l-/i)] 



hi ^ Sjl-h) 
ho ho 



boundary layer taken in the region before the jump, and 



+ Sul (/i - /a) = - (hi - hi) (8) . 

2 as such, its value must be less than hi . This inconsistency 

As we have done for the inviscid case, we use the same def- is indication enough that e is not zero, and therefore the 

inition of H and JF, and obtain an expression from Eq.(|Hl transition is not continuous. 
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Fig. 3. Plot of ${e) versus e. The roots of e as given by Eg. llUII . 

are to be found for ^(e) = 0. Both plots have been drawn 
for S/hi — 0.5, with the continuous curve representing the 
parabolic profile, and the dotted curve representing the profile 
for the Couette flow. 

For a more methodical evaluation of the roots of e, 
we will have to solve Ea. (|l()(l . To that end, for notational 
convenience we write the third degree expression on the 
left-hand side of Ea. Hll)|) as and then solve for <?(e) = 
0. To find the roots of this cubic equation, it would be 
necessary to eliminate the second degree term in e, by 
the suitable substitution e = C — 1. This will then render 
Eg. pUfl in the standard form as C'^ + T^C +2 = 0, where 
P = -3 + 2(1 - h'jS/hi and Q = 2[1 - (1 - l2)5/hi + 
(1 - Ii)^{S/hi)^]. The discriminant, V, is given by P = 
(QV4) + {VV27). 

We choose 0.5 as a fiduciary value for 6/ hi, and using 
this number for both the parabolic profile and the profile 
for the Couette flow, we find V < 0. This must then im- 
ply that there should be three real roots of e for — 0. 
This is very much in keeping with the fact that the func- 



tion has two real turning points, which can be ob- 
tained from the condition 'P'{e) = 0. These points are at 
e = -1 ± v^l - (2/3)(l - l2)S/hi. Further, when e = 0, 
the function has a negative value. This, in conjunc- 
tion with the fact that ^(e) — > +oo when e — > +oo, 
could only mean that at most there should be only one 
positive real root of e, a conclusion that has been clearly 
illustrated in Fig|21 in which we have plotted ^(e) against 
e. For both the Couette and the parabolic profiles, three 
roots are to be found for 'P{e) — 0, and in each case only 
one of these roots is positive. A numerical evaluation by 
the bisection method shows that for a parabolic profile the 
positive root of e is about 0.15, while for the Couette flow 
it is 0.17. For both these cases we have also assured our- 
selves that even for the limiting case of 6 /hi being unity, e 
will have a positive root. The inescapable conclusion there- 
fore is that the jump is a first-order transition, indepen- 
dent of the vertical velocity profile, because at the critical 
Froude number (J- — 1) the jump has a finite non-zero 
value. The present analysis drives home the point that 
when dissipation is included in the conventional control 
volume analysis, the second-order transition immediately 
changes to a first-order transition. 

3 A Hydrodynamic Analysis from the Steady 
Flow Equations 

We consider an open rectangular channel of width L (as 
shown schematically in Fig^ in which a liquid is flowing 
in streamline motion. An arrangement is made such that 
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the liquid flows down an inclined channel and starts its gible. In the boundary-layer approximation it is assumed 

one dimensional motion in the horizontal plane from the that the vertical velocity w{x, z) is very small compared 

edge AB. The x axis is chosen in the direction of flow, the with the horizontal component u{x, z). Furthermore, the 

y axis along the width of the channel and the z axis in the variation of u{x, z) is much faster in the z direction as 

vertical upward direction. At any point, the velocity of the compared with the x direction. 

liquid is, in general, a function of the coordinates x, y and At this stage we follow Bohr et al. 9 to do an averag- 

z. In our present arrangement the width of the channel ing of the flow variables over the z direction. We define 

is sufficiently large (about 9.1cm) when compared with rH^) 

(^(^'^)) = / '<P{x,z)dz (14) 

the height of the liquid layer (which is of the order of a 



Jo 



few millimetres), and we can therefore assume that there where the averaged quantity in the angled brackets be- 

is no variation along the y axis for the velocity. On the comes a function of x only. Under the assumption of a flat 

other hand, both its x component, u, and z component, w, free surface such that ■w{x, z) = at z = h{x), and along 

would have a spatial dependence in the form u = u{x, z) with the use of Eq.(ITTl), we can readily show that 
and w = w(2;,z). / du\ / du\ 



dz I \ Qx 

For an incompressible fluid, the local continuity equa- 
tion gives, ^^'^ 

du dw _ /^\ = 

a^+aJ-O (11) \dz2/ 

and for a not very viscous liquid (e.g. water) the Navier- We now make the approximations. 



(16) 



z=0 



Stokes equation in the boundary-layer approximation USI / du\ d(u) 

)=«M4r (17) 

gives 



dx / dx 



du ^ du dh ^ d'^u 

dx dz dx dz'^ 

where h = h(x) is the height of the liquid layer at a 



and 

du 
'd~z 



z=0 



distance x, and v is the kinematic viscosity. The bound- where a and (3 are numbers of 0(1) and they depend upon 

ary conditions of the flow are w(x,0) = w(x,0) = 0, and the velocity proflle. For a parabohc profile a = 6/5 and 

du/dz = at z = h{x). In addition to this, the condition /? = 3. It is easy to check that the parameters a and (3 are 

for constant volume fiux gives strictly constants when we make, following Watson !B|, the 

i,h(x) reasonable scaling ansatz that m(x, z) = U{x)F[z/h{x)\. 
L I u{x,z)dz = Q (13) 

•^0 The fact that our experiment will determine a combina- 

We have assumed here that the shearing stress is zero at tion of a and (3 is important regarding the nature of the 

the free surface z — h{x), since the viscosity of air is negli- vertical profile of the velocity field over which we average. 
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With the above approximations and identity, and writ- where C is a constant of integration. The position Xj of the 
ing (u) as v, Ea.l(T^ becomes jump is obtained by setting h ^ hj m Eq.|(5lJ- This yields 

2„„i!: = -,ii_^„" ,19) ' ii\4j [l) 

dx ' ax h'' 

where C is some other constant. 
From Eq.lHSJ we get Lvh = Q, which we can use to elim- The complete profile of h{x) is described by Eq.^, 
inate v from Ea. lfT^ and get but it should be noted that it gives multiple h{x) for a 

given X. However, for values of h{x) which are even mod- 

dh _ Q 

~ y^^) erately less that /ij, the profile according to Ea. pifl can 

be approximated as 



2 



h{x)^['^^\x (23) 



The derivation of the above equation may look somewhat 

heuristic, but it has been much more systematic to the \2a Q 

extent that it avoids the inclusion of an artificial friction ^hich tells us that for small x the height of the liquid 

loss. The structure of this equation, however, is quite sim- layer increases linearly. This feature is markedly different 

ilar to what is well known in hydraulic jump literature ;6,, ^om the case of the radial spread of a liquid stream on 

to derive which, the conventional recourse has been to use a horizontal plate in which the height gradually decreases 

the Bernoulh equation supplemented by a friction loss. In t^ie region within the jump [T]. 

our derivation ofEq.®, the approximations made on the With the help of Eq.jai) we are now in a position 

Navier-Stokes equation have been clearly delineated, and to get a physical picture of the jump. We introduce the 

the resulting profile has been seen to be sensitive to the ve- dimensionless variables H = h/h^ and X^x/i-^ in which, 

locity field in a fashion that can be experimentally probed. ^ J =^i~C- We thus obtain a dimensionless form of Eq. ^ 

This equation shows that li v = {), then h{x) would be a as 

constant. Therefore, without viscosity the hydraulic jump ^ ^ = ^{X — D) (24) 

can be obtained only if we explicitly seek a solution with with D being a dimensionless constant. Differentiating 

the jump extraneously imposed, which is the way it has Eq.(|23I) will readily show that at = 1, the function 

been conventionally treated in text books. The jump oc- H{X) will have a vertical tangent — something that can 

curs when dh/dx displays singular behaviour at /i = /ij be conceived of as a discontinuity in the physical flow. Our 

such that = 2a{Q/LY. The advantage of Eq.(|2n|) is analysis is restricted to the range X > Q. The point along 

that it can be exactly integrated. This gives, the X-axis, where the flow will encounter this discontinu- 
ity, will be determined by the inner boundary condition. 

9^ [^3 ~ 1" j ^ ^'^J^^ + (21) To show this, we set the condition H = at an arbitrary 
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Fig. 4. Fitting the inner and outer solutions via a shock. The 
inner solution corresponds to X < 1, while the outer solution 
corresponds to X > 1. Each solution has been determined by 
a different boundary condition. The vertical dotted line (the 
shock) joins the two solutions. 

X = X[n, to get D = X[^. When X = 1+Xin, we will have 
H = 1, which gives us a clear indication that the value of 
will determine the position of the discontinuity. 

A generic inner boundary condition is that _ff = at 
A" = (ignoring a small initial non-zero value of -ff at AT = 
0), giving us A"in = D ^ 0. The solution corresponding to 
this particular boundary condition is represented by the 
lower branch in Fig^] In the dimensionless notation that 
we have introduced, this means a profile proceeding from 
X = to X = 1 and rising from H = to H = 1. Since 
the inner solution diverges infinitely at AT = 1, the flow 
has to physically go beyond this point by the fitting of 
the lower branch to the upper branch in Fig^ through 
the discontinuity. We treat this discontinuity as the jump 
in the flow. 



The integration constant D for the upper branch is 
to be fixed by the outer boundary condition. This is de- 
termined by the physical requirement that H = 1 (in- 
finite slope) at A' = ATc, with Xc corresponding to a 
position slightly left to the edge of the channel where 
the liquid falls off. We cannot extend the solution to the 
edge due to singular fiow at the outlet. For the outer 
boundary condition stated above we will get the solution 
4H - ^ 3 {X - Xc + 1). At the position of the jump, 
i.e. at X = 1, the outer solution will give the magnitude 
of the jump, H], making it evident that Hj has a depen- 
dence on Xc- It should be a worthwhile exercise to study 
this dependence. 

The roots of Hj can be determined by solving the bi- 
quadratic equation given by Hj — 4Hj + 3 (2 — A'c) = 0. To 
that end, we will first have to carry out a transformation 
with the help of a new variable 77, to write (iJj + 77)^ = 
2f]H] +if + AH} - 3 (2 - Xc), and then require that by 
a suitable choice of 77, the right hand side will be ren- 
dered a perfect square. This requirement will mean that 
the discriminant of the quadratic in iJj on the right hand 
side should vanish, to ultimately yield the auxiliary cubic 
equation if - 877 (2 - Xc) - 2 = 0. 

The discriminant, A, of this cubic equation will be 
given by Z\ = 1 — (2 — X^f ^ and it is quite easy to see 
that for Xc > 1, there will always be a positive value 
for A. Therefore, 77 can have only one real root, 7yo, given 
by ?7o = 1 + \/^ + 1 - ^) ■ To solve for i/j. 



we now use this value of 770 in the biquadratic equation 
(i/j + 770) = ±^2770 + Of course, there wiU be 
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four roots of Hj, but the real roots should correspond 4 Time-dependence in the Channel Flow 

to the choice of the positive sign. Solving the relevant 
quadratic equation and choosing to keep only the phys- 
ically meaningful positive root, will give us the solution of 
iJj as 



So far our analysis has been carried out in terms of the 
steady flow equations only. We now consider explicit time- 
dependence in both the equation of continuity and the 
Navier-Stokes equation. This will necessitate the time- 
2 V V^o/ 2 ^ ^ dependent generalisation of both the flow variables /i, the 

An interesting conclusion that can be derived from the local height of the flow, and v, the vertically integrated av- 
relation of the jump height above is that for Xc ^ I, the erage flow velocity. The resultant dynamic equations may 
maximum height of the jump will asymptotically be given then be written as 
byffJ^(3Xe)^/^ dh d 

As a specific case we set Xc = 2, and then determining 
the values of both A and r^o, we see that i7j = 4^/'^. This 

dv ^ dv ^ dh V 
result could alternatively be derived directly from Eq. H24II dt dx dx 

by the use of the boundary condition, H = I, aX X = 2. respectively. 

This will give D = 1, and for this particular boundary For the purpose of carrying out a linear stability anal- 
condition, the resulting outer solution has been plotted in ysis of the flow in real time, it will be convenient for us 
FigH In this instance the fractional change in height at to work with a new variable which we define as / = vh. 
the jump position would be 4^/^ - 1 ~ 0.59. The new variable / can be physically associated with the 
The fitting of the upper and the lower branches (each time-dependent volumetric flow rate, and its steady solu- 
determined by a different boundary condition) has to be tion, as can be seen from Eq.(|lll), is a constant. We use 
done via a shock, by requiring that the mass and momen- solutions of the form v{x, t) — vq{x)+v {x, t) and h{x, t) — 
tum fluxes remain unchanged through the discontinuity — h^{x) + h' {x,t), in which the subscript indicates steady 
a condition that may be derived from the inviscid equa- solutions, while the primed quantities are time-dependent 
tions. From Eqs. (HJ and @ , with both hi and /12 scaled by perturbations about those steady solutions. Linearising in 
/ij to give Hi and H2, we get what Bohr et al. call the these fluctuating quantities gives us the fluctuation of / 
jumphne H1H2 {Hi + H2) = a'\ The jump takes place about its constant steady value /q = voho = Q/L, as 



when H of the viscous solution in the lower branch of Fig2] 
becomes equal to Hi of the inviscid jumpline for a given 
value of a. In general, the shock fitting leads to a jump 



/' = voh' + hov' (28) 



In terms of /', we can then write from Ea. (|26|l . 



dh' _ df 

position actually slightly left of X = I. -q^ — ~~q^ (^^) 
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and this in turn, along with Ea. H28|l . gives, 



same wave equation as for a scalar field in curved space- 
time 14 . For an inviscid, incompressible and irrotational 
flow, these authors have prescribed the flow velocity to be 
A further partial derivative of Eq. m with respect to time g^^Hent of a scalar potential. Perturbing this potential 



dv' _ 1 
dt Hq 



dt 



vo fdr_ 

dx 



(30) 



yields 





d 






d 








dt 


ho 




"-at 


ho 





(31) 

The significance of the form in which we have kept Eq. H31|) , 
will be apparent soon. Linearising in the perturbed quan- 
tities in Ea. (|77|) gives 



dv 



d 



dh' 



h' 



-d^ + Q-,M+9^ = -^(v'-2voj^^ 

which in turn, upon partially differentiating with respect 
to t, yields 



2„/ 



d 



dt'^ dx 



Wo 



d_i/_ 
dt 



d fdh' 
t 

dv' 



^ dx V dt 



2^ ( — 

dt ho V dt 



In Ea. H33|l above, we substitute for the first and the second 
order time derivatives of h' and v' from Eas. H29|l to (|31|) . 
This will lead to the result 



d 






d 




m] 


d 


Vo 


dt 


ho 




^dt 


ho 




dx 


ho 



d_ 

dx 



Vo - '^'^ ^ 



/in \ dt "° dx 



about its background value, under the restricted condition 
of the flow height being a constant, leads to an effective 
metric of the flow, in which the velocity of gravity waves 
replace the speed of sound in sonic analogs that faithfully 
reflect features seen in general relativistic studies |14II15| . 

We now proceed to demonstrate that our perturbative 

(32) analysis of what is essentially a dissipative system (since 
it includes viscosity), will, in its inviscid limit, deliver the 
same metric obtained by Schiitzhold and Unruh from their 
purely inviscid model. It must be stressed here that our 
choice of perturbing the flow rate is paticularly expedient, 

(33) since conservation of matter holds good even in a sys- 
tem that undergoes viscous dissipation. It is to be further 
noted that the background velocity and flow height in our 
treatment are in general functions of space and not con- 
stants. Having made these observations to emphasise the 
greater generality and exactitude of our approach, we can 

(34) then extract the inviscid terms from Ea. (|34|l by setting 
= 0, to ultimately render these terms into a compact 



formulation that looks like 

(f^'d^f) = 



(35) 



At this juncture it should be most instructive for us to 
examine Ea. (|34|) in its inviscid limit, i.e. when v = 0. In 
connection with this, it should also be very much worth 
our while to consider some recent studies jl4lll5lll6) which 

have revealed a close analogy between the propagation of a in which, we make the Greek indices run from to 1, with 
wave in a moving fluid and of light in curved space-time. In the idcntiflcation that stands for t and 1 stands for x. An 
particular Schiitzhold and Unruh have shown how gravity inspection of the terms in the left hand side of Eq. (|34|l — 
waves in a shallow layer of liquid are governed by the all of them divided by the constant g — will then enable 
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us to construct the symmetric matrix governed by inviscid conditions, and that the basic prop- 

1 / 1 vo \ erties of gravity waves wiU not be drasticaUy affected. In 

^'"-zri (36) 

^ " \ Wo Vq — gho / our presentation of the perturbative analysis, we have sys- 

Now in terms of the metric g^t^, the d'Alembertian for a tematically included viscosity in our governing equations, 
scalar in curved space, is given by JS] which has finally led to Ea. H34|l . It is quite evident that 

1 the presence of viscosity disrupts the precise symmetry of 

A7A^-=a^(^/^g^''(}.V) (37) 

V S the inviscid conditions described Ea. (|35|l . This obviously 

in which g^'^ is the inverse of the matrix implied by g^^i/. implies that the clear-cut horizon condition that one ob- 

Under the equivalence that f^"^ — ^—gg'^", and therefore, tains from the inviscid limit, will be affected. However, 

g = det (f^'^), we can immediately set down an effective this effect, for small viscosity, as Schiitzhold and Unruh 

metric have argued, cannot be too drastic. This is very much in 

— (3g) conformity with our observation in the previous section 

\vo vi~gho/ 

that fitting the inner and outer solutions through a shock 

that is entirely identical to the one obtained by Schiitzhold 

will shift the position of the jump only slightly. One way or 

and Unruh, following whom, the inverse effective metric, 

the other, the most important feature to emerge from the 

g^^i can also be easily obtained from Eq.JJHJ- This shall 

analogy of white hole horizon shall remain qualitatively 

identify Vq — ghQ as the ergosphere condition on the hori- 

unchanged, namely, that a disturbance propagating up- 

zon of either a black hole or a white hole, depending on the 

stream from the subcritical flow region (where Vq < gho) 

direction of the flow. In the case of the two-dimensional 

cannot penetrate through the horizon (where Wq — ghg) 

circular hydraulic jump, as Volovik has pointed out |16| . 

into the supercritical region of the flow (where Vq > gho). 

the jump condition can be closely related to the horizon of 

As we shall see shortly, this property of the flow will have 

a white hole, a surface that nothing can penetrate. This 
analogy is entirely apt for the channel flow that we are 
studying here, considering the direction in which the flow 
proceeds. 

-r , . „. ...in 1 rr,i.i the bchaviour of the perturbation f. To that end we go 
Introduction of viscosity m the now does attect the ide- 

back to Ea. (|34|l and recast it in a slightly altered form 
that looks like 



a very crucial bearing on a physical picture that we shall 
construct to explain the formation of the hydraulic jump. 

For our purposes it should also be important to study 



alised inviscid conditions in the vicinity of the white hole 
horizon. Schiitzhold and Unruh themselves have treated 



viscosity as a small adjunct effect on the inviscid flow, q2 j^i d f df'\ Id 
and have concluded that although viscosity will introduce 



dt'^ dx \ dt I Vq dx 



df 

vo (v^o - gho) -g^ 



V (df df'\ 

a thin boundary layer, the flow outside it shall very well be ~ \ dt ^ '^^^~dx J ^"^^^ 
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Using a solution of the type f'{x,t) = p{x) exp{—iojt) in Ea. (|40|) leads to 

Eq . 1391) ffives an expression that is to be further multiplied p 2 r / \ j 

2 f 2, , ■ f voP , o f I dp 
u! / vqp ax + \ljjv I —p^ ax — iv I ( — ) p— ax 

by vqp throughout, to finally deliver a quadratic equation •' 



2 

. I ,2 , \ I 

in Cij that is of the form 



/iq J \hQ J dx 

vo{vl-gho) (^^) dx = (41) 



cj^ (vqp'^) 



— (VaP) + v—^ 



dx 



which is a result that has been arrived at by carrying 



d 

ax 



y2p / (^p\ out the integration by parts, and then by imposing the 

hn \dx. 



3-^(^1=0(40) 



requirement that all the integrated "surface" terms vanish 

To have any idea of how the perturbation behaves in time, a-t the two boundaries. 

we treat it as a standing wave. For that purpose we will Under inviscid conditions, lu will have a purely real 
have to integrate the above equation between two chosen solution, and therefore the perturbation will be a stand- 
boundaries, at which the perturbation will be constrained ing wave with a constant amplitude in time. However, the 
to vanish at all times. Between these two boundaries the dissipative presence of viscosity will result in the pertur- 
flow should be continuous. Since we are aware that the bation being damped out, and will restore the system to- 
jump itself is a discontinuity in the flow, we will have to wards a stable configuration. The consequent time-decay 
choose the boundaries to be on one side of the jump only, of the amplitude of the perturbation will be of the form 
although Eq. (gni itself holds true over the entire range of exp{~-vt/2hl). This also gives a time scale for viscous 
the flow. We have already acquainted ourselves with the dissipation, which, to an order-of-magnitude, is given by 
fact that a perturbation in the subcritical region will re- tvisc ~ ^o/-- 

main confined to this region only. Besides, in this region It is now important to appreciate that viscous drag 
the flow would have entirely lost its laminar character, in the fluid will also dissipatively slow down the flow on 
and therefore, would be most suited for us to derive some this very same time scale ivisc- The information of an ad- 
physical insight about the behaviour of the perturbation vanced layer of fluid slowing down has to propagate up- 
and the influence of viscosity on that. Therefore, we con- stream to preserve the smooth continuity of the fluid flow, 
fine our analysis to the subcritical region of the fiow only. This propagation, however, cannot happen any faster than 
As for the boundaries of the perturbation, one of them the speed of the surface gravity waves, {ghoY^^, and in 
can be the outer boundary of the steady flow itself, while the region where vq > (ghoY^"^, i.e. in the supercritical 
the inner boundary may be chosen to be very close to the region, no information therefore can propagate upstream 
jump. In this regard we treat the jump itself as a bound- |11| . This is also what we should entirely expect from the 
ary wall where all velocity and height fluctuations decay modelling of the jump as the horizon of an impenetrable 
out completely. Under these conditions an integration of white hole. So a stream of fluid that has arrived later, af- 
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ter having passed through the supercritical region, moves the fluid, and what we see is a sudden increase in the 
on ahead, yet unhindered and uninformed, till its speed height of the fluid layer — a hydrauhc jump. This gives a 
becomes comparable with the speed of the surface gravity conceivable physical basis for understanding how crucial a 
waves, and only then does any knowledge about a "bar- factor viscosity is behind the formation of the jump, and 
rier" ahead catches up with the fluid. We now define a in connection with this physical picture, it is also worth- 
dynamic time scale, idyn ^ x/vq, on which the bulk flow while to ponder the possibility that the viscosity-induced 
proceeds. If we set ivisc — ^dyn with the additional require- boundary layer of the flow gradually increases in thickness 
ment vq — (g/io)^^^, we get the condition for the "news" and the jump occurs at that distance, where the bound- 
of the viscous slowing down finally catching up with the ary layer pervades the entire height of the thin layer of 
bulk flow itself. With the additional constraint that we the flowing fluid, i.e. from z = to z = h{x). 

have from the continuity equation, i.e. vohn = QIL, we ttti -i j n- j.i • r j.i i • i • j. 

•' ^ ' 'J " -K/ I While dweUmg on the issue of the physical picture 

can then derive a scaling relation for length which is en- c ■ c ^- ^ ij^jj 

° oi the jump formation, we are tempted to adduce a re- 

tirely identical to the scaling dependence of the position i.j, i-i i ^- i- 

t, f f lated astrophysical analogy : the formation of hot spots m 

of the hydraulic jump, ^ ■ • i ^- t ^ ■ i 

J galaxies. Gaseous jets emanating from galaxies encounter 

5/3 



— J V (42) resistence from the intergalactic medium. As a result the 

1 J?™ li _Li i 1 111 tip of the iet is slowed down in comparison with its bulk, 

as given by £jq.(|22|) — a result that we have already ob- f ■> 

, . J r i J r i £1 £• This causes energy to accumulate at the tip, and this is a 

tamed from our study of the stationary flow equations. ^■^ 

mi r i-1 i i-1 i- r i likely explanation for the formation of a hot spot. The gas 

1 he crux of the argument that emerges Irom our anal- j f t- o 

• • j-u J- r j-u r -t- r j-T, T, J r • -ti Aow in the jet proceeds at supersonic speeds, but on com- 

ysis IS that lor the lormation ol the hydraulic jump, the j ^ j- t- j 

j^,. , ^1 1 ing close to the hot spot, the flow experiences an abrupt 

two time scales, ivisc and idyn, would have to match each ° 

1 1 1 ii 1 1 -r- • 1 i -i deceleration. The information of this sudden braking is 

other closely when the l^roude number, y- , is close to unity. ° 

TT J ii i-i- 1 rfl -J • ■ 1 i- 1 • not conveyed upstream, since no sound wave can move 

Under these conditions, a layer of fluid arriving lately is j f i 

, 11 1 rn-1 • against the supersonic flow. This causes a shock wave to 

conlrontcd with a barrier lormed by a layer ol nuid moving ° 

, 1 1 i 1 mi • 1 1 -1 form across the iet 1171 . In this particular astrophysical 
ahead with an abrupt slowness. 1 his slowly moving layer 

rn-jn j , T ■,• j-j- plcturc, ouc might discern much similarity of principle 

ol fluid flowed past earlier m time, and at lar distances ^ ' <= j t- t- 

., 1 , i 1 J -J ui 1^ • J o- with our physical arguments behind the formation of the 

it has been retarded considerably by viscous drag, bmce ° 

. ii . ii 11 ,1 . , „ hydraulic iump in the channel flow, 

m this situation there could not be an indefinite accumu- o l 

lation of the fluid, and since continuity of the fluid flow Another significant point of which we should like to 

has to be preserved, the newly arrived fluid layer slides make mention is that instead of Ea. (|27|) we could have 

over the earlier viscosity dragged slowly moving layer of chosen to use the dynamic generalisation of Eq. (|19|l with 
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the dimensionless constants a and /? included, but to de- small extraneous addition to the flow should not too dras- 
rive the particular form of the perturbation equation as tically affect our observations. In any case, preserving the 
given by Ea. (|34|l . we would have to set a — 1/2. This ar- laminarity of the flow should be necessary for the unam- 
gument possibly has an important bearing on the issue of biguous recording of data. For that purpose, the height of 
the velocity profile of the flow. the water layer has been measured with a travelling micro- 

scope. The flow has been viewed through the transparent 
5 Experimental Results perspec wall of the horizontal channel. At each value of 

a;, we noted down the vertical positions of the bottom of 



A relatively simple experiment, using water, was carried 
out to substantiate our theoretical propositions. Our ob- 
jectives were two-fold. First, to verify the linear growth of 
the surface height of the flow at small distances, and sec- 
ondly, to verify the scaling of the jump position as given 
by Ea. H42|) . It has been satisfying to note that our theory 
— presented heuristically in parts — has very much been 
borne out within the limited objectives of our experiment. 

The experimental apparatus has already been illus- 
trated schematically in Fig^ The water flows in a stream- 
line motion through an open rectangular channel of width 
L. An arrangement has been made such that the water will 
flow down an inclined channel, with an inclination of 60°, To estimate the volumetric flow rate, we adopted the 



the flow in contact with the bed of the channel, and of 
the free surface of the flowing water (both of which were 
easily identiflable in the fleld of view of the microscope). 
The difference of the two readings gave h(x). 

Viewed from the top of the channel, the jump itself 
has been seen to present a curved front across the width 
(i.e. along the y axis) of the channel. This is because the 
flow in contact with the boundary walls of the channel (at 
y — Q and y — L respectively) has been dragged down 
by viscosity. The position of the jump, xj, is actually an 
average value measured over this lateral curved proflle. 



and then from the edge AB, will start its one-dimensional simple recourse of collecting the water falling off at the 

motion in the horizontal channel whose length is 70 cm o^^r edge of the channel, and then of measuring the vol- 

and width is 9.1 cm. A rectangular box is attached to the of the water collected for various intervals of time, 

base of the inclined channel. Water comes from a tap to The average of all these readings has been taken to deter- 

the rectangular box. There is a sHt in the box through ^ine Q. The steadiness of the flow has also been conflrmed 

which water flows down the channel. The purpose of the by this approach. Values of X,, at the outer edge of the 

inclined channel is to make the flow laminar from the be- channel (discussed in Section EJ, range between 1.5 and 

ginning of the horizontal channel. Although this will in- 2.8 in our experimental set-up. 

troduce a small horizontal component to the flow, in so far We observed a very slow rise in h{x) for a while and 

as we would be interested in the scaling relations only, any then a major jump. That the rising proflle of h{x) will be 
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Fig. 5. The height profile h{x) before the jump for different 
values of Q, the volumetric flow rate. 
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Fig. 6. The position Xj of the jump for different volumetric 
flow rates Q. 

linear, as we might rightly expect on the basis of Ea. H23() . 
has been shown in Fig|31for two different values of Q. We 
also show the scaling dependence of the jump position Xj 
on Q^^^ in FigEl The predicted linearity from Ea. (|42ll 
has been depicted quite clearly here. It is obvious that 
once the 5/3 power dependence of xj on Q has been es- 
tablished experimentally, the dependence on and g, as 
Eg. 1)42(1 shows, must follow even on the basis of elemen- 




Fig. 7. A sideview of jump region with the flow proceeding 
from the left to the right. The flow appears black in the photo- 
graph, because it has been coloured with a red dye. The jump 
is clearly discernable, as also are the growth and decay profiles 
of the surface height, before and after the jump, respectively. 

tary dimensional considerations. So our theoretical scaling 
law stands vindicated by our simple experiment. 

More to this point, we also furnish two photographs of 
the cross-section of the flow with the jump included. We 
show a long distance snapshot of the flow in Figd The 
flow has been dyed red to make it more prominent, when 
viewed through the transparent perspec wall. The jump 
is very much discernable in this photograph, but we must 
also draw attention to the slow linearised growth of the 
height of the flow much before the jump, followed by its 
much more rapid growth immediately in front of the jump. 
This is further followed by a small decrease in the flow level 
in the region beyond the jump. Qualitatively this is what 
we should expect on the basis of the plot in Fig0] That 
the variation of the flow height we see in the photograph, 
is not as pronounced as Fig^ would impress upon us, is 



16 



Subhendu B. Singha, Jayanta K. Bhattacharjee, Arnab K. Ray: Hydraulic Jump in One-dimensional Flow 




Q 




^2 


£ = H - 1 


(cm^ s^-^) 


(cm) 


(cm) 




66 


0.415 


0.590 


0.42 


91 


0.430 


0.619 


0.44 


121 


0.445 


0.648 


0.46 


148 


0.573 


0.856 


0.49 


152 


0.590 


0.890 


0.51 


162 


0.621 


0.940 


0.51 



Fig. 8. A close-up sideview of the jump region. 

because of the fact that the axes in Fig0|have been scaled 
in terms of the flow constants, while the variation of the 
flow height in the actual photograph proceeds on the scale 
of CGS units. The jump appears much more prominently 
in the close-up view of the jump region, shown in Fig|Sl 

From Tabled it is seen that the jump remains almost 
constant within an experimental error of about 6%. The 
combined effect of the theoretical analysis and the exper- 
imental results leads us to believe that the magnitude of 
the jump in a shallow layer flow will show no drastic vari- 
ation. 

The z profile of the velocity field u[x,z) also calls for 
some comments. The slope in Ea. (|22ll depends on the 
numbers a. and /3 according to the combination a^^^ / [3, 
which, for a parabolic profile (a — 6/5 and /3 = 3) is 
about 0.425. This ultimately gives a theoretical estimate 
of the slope to be approximately 0.204. But the experi- 
mental data as plotted in FigElshow that the actual slope 
(nearly 0.007) is much smaller. This leads us to infer that 



Table 1. The magnitude of the jump e, at different volumetric 
flow rates Q. 

the profile is much steeper than parabolic near the plate 
and much flatter near the free surface. 

In the experiment we have also seen that the flow be- 
comes turbulent after the jump, with the formation of vor- 
tices. We have not made any measurement for this region 
as we have been interested in the laminar flow only. How- 
ever, we have noticed that for low volumetric flow rates, 
turbulence in the subcritical flow region dies down ap- 
preciably, as the flow proceeds downstream. On the other 
hand, for high flow rates, turbulence is seen to be sus- 
tained right upto the outer boundary of the flow. Quali- 
tatively this is what it should be. Turbulent fluctuations 
derive their energy from the mean flow. If the mean flow is 
more energetic, such as it should be for higher flow rates, 
then turbulent effects will linger in the flow that much 
longer. Regarding this issue, it should be possible for us, 
at least in an order-of-magnitude sense, to make a theo- 
retical estimate of the Reynolds number of the turbulent 
flow, immediately after the jump. In the shallow layer flow, 
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the largest possible turbulent eddies should have a char- 
acteristic length scale that should at most be of the order 
of the flow height h, in the subcritical region. The char- 
acteristic turnover velocity of the eddies should likewise 
be of the order of (but less than) the velocity of surface 
gravity waves, (g/i)^/^. The Reynolds number of the flow, 
TZc, should therefore be given by 7^c - v-^igh^fl'^. In 
our experiment, typical values of the flow constants g and 
V would be 1000 cm sec~^ and 10^^ cm^ sec~^, respec- 
tively. From Tabled for various values of Q, we may have 
an estimate of the characteristic values of the flow height h 
immediately after the jump. These measurements should 
then typically give TZc ^ 1000, an estimate, whose direct 
verification, however, would be beyond the scope of our 
experiment. 
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